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Recent work on gravitational geons is extended to examine the stability 
properties of gravitational and electromagnetic geon constructs. All types 
of geons must possess the property of regularity, self-consistency and quasi- 
stability on a time-scale much longer than the period of the comprising 
waves. Standard perturbation theory, modified to accommodate time- 
averaged fields, is used to test the requirement of quasi-stability. It is found 
that the modified perturbation theory results in an internal inconsistency. 
The time-scale of evolution is found to be of the same order in magnitude 
as the period of the comprising waves. This contradicts the requirement 
of slow evolution. Thus not all of the requirements for the existence of 
electromagnetic or gravitational geons are met though perturbation theory. 
From this result it cannot be concluded that an electromagnetic or a grav- 
itational geon is a viable entity. The broader implications of the result are 
discussed with particular reference to the problem of gravitational energy. 
PACS 04.25.-g 



1 



I. Introduction 

The examination of the basic properties of the gravitational field as compared to other physical 
fields has concentrated around the recently revived study of gravitational gcons. The concept of 
a structure comprised of electromagnetic waves held together by its own gravitational attraction 
was first conceived by Wheeler. 1 The extension of this idea using only gravitational waves was first 
studied by Brill and Hartle. 2 Their approach was to consider a strongly curved static 'background 
geometry' 7^ on top of which a small ripple h^v resided, satisfying a linear wave equation. The wave 
frequency was assumed to be so high as to create a sufficiently large effective energy density which 
served as the source of the background 7^, taken to be spherically symmetric on a time-average. 
They claimed to have found a solution with a flat-space spherical interior, a Schwarzschild exterior 
and a thin shell separation meant to be created by high-frequency gravitational waves. With the 
mass M identified from the exterior metric, there would follow an unambiguous realization of the 
gravitational geon as described above. It has since been shown 3-5 that the Brill and Hartle model 
docs not implement the properties of high-frequency waves, nor can the space-time be taken as 
singularity- free. 

It was proposed by Cooperstock, Faraoni and Perry 3-5 (henceforth referred to as CFP) that a 
satisfactory gravitational geon model must be constructed and solved in a manner similar to that 
of Wheeler's 1 electromagnetic geon model. Such a model necessarily requires firstly that the Ein- 
stein field equations be solved in a self-consistent manner while satisfying the regularity conditions. 
Secondly, it is required that the configuration represented by the metric 7^ be quasi-stable over 
a time-scale much larger than the typical period of its gravitational wave constituents (i.e. a geon 
must maintain its bounded configuration for a sufficient length of time, for otherwise it would not 
be possible to attribute a structural form to the gravitational geon). Thirdly, it is required that 
the gravitational field becomes asymptotically flat at spatial infinity. Thus a gravitational ( electro- 
magnetic ) geon is a bounded configuration of gravitational (electromagnetic) waves whose gravity is 
sufficiently strong to keep them confined on a time-scale long compared to the characteristic com- 
posing wave period. For the gravitational case, it is required that no matter or fields other than the 
gravitational field be present. 

Through a series of papers, it was established by Anderson and Brill 6 and by CFP that in the 
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high-frequency approximation for a static background metric, the gravitational geon problem and 
the electromagnetic geon problem are governed by the same set of ordinary differential equations 
(ODEs) and boundary conditions. These equations are satisfactory for considering the regularity 
and self-consistency aspects of the geon problem but not the evolution in time. Any solutions to 
these equations are necessarily equilibrium solutions since the background metric is assumed static. 
Admissible equilibrium solutions satisfying the boundary conditions have been shown to exist. 6 
This paper provides an expanded study of the gravitational and electromagnetic geon problem 
with particular emphasis upon the dynamic evolution of geon constructs. Section II re-examines the 
numerical solutions to the equations studied in Ref. 1 and extends the analysis to obtain information 
on the "stability" of these solutions with respect to perturbations of the amplitude eigenvalues. The 
word "stability" used in Sec. II should be viewed in the context of the boundary conditions of a 
spatial variable, not a dynamic time variable. Both the analytic behaviour at spatial infinity and the 
aforementioned stability properties are found by constructing a phase portrait of the ODEs. It is 
explicitly shown that only unstable equilibrium solutions are possible with respect to perturbations 
of the amplitude eigenvalues. This result serves to suggest further study of the dynamics (time- 
evolution) of geon constructs. 

The evolution in time of electromagnetic geon solutions is studied in Sec. III. The evolution of 
the electromagnetic geon is studied instead of the gravitational geon because of the relative ease in 
computation for the former. With sufficiently high-frequency electromagnetic waves, the results are 
expected to apply equally well to the gravitational case. The method used is standard perturbation 
theory modified to accommodate time-averaged fields. The evolution is achieved by applying a small 
amplitude time-dependent perturbation to an equilibrium solution and simultaneously solving for the 
time-dependence of the background metric functions. The problem of time-averaging the background 
metric functions can only be done in a meaningful way if it is assumed that the characteristic 
frequency of the perturbations vary on a time-scale much longer than that of the waves comprising 
the electromagnetic geon. This is in accordance with the requirement that the background metric 
be quasi-stable. However, the results of the analysis show that the perturbations must vary on 
the same time-scale as the constituent waves. This is a contradiction to the original assumption. 
Hence an internal inconsistency exists when applying perturbation theory to the geon problem. In 
Sec. IV, the possible interpretations and ramifications of this result are discussed. Since not all of the 
requirements for existence of a geon are met, it is not possible to conclude that an electromagnetic 
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gcon or a gravitational geon is a viable construct. The conclusions are presented in Sec. V. 

II. Phase Space Analysis 

In the high-frequency approximation, the gravitational and electromagnetic geon problem for a 
static background metric (on time-average) reduces to the same set of ordinary differential equations 
(ODEs) and boundary conditions given by 3-6 

4>" + jk<t> = Q, (1) 
k< = -cf> 2 , (2) 
.?' = 3-fc- 2 (l + ^ 2 ), (3) 

and 

4>(x) — > , k(x) — > 1 and j{x) — > — oo for x — > — oo, 

(4) 

4>{x) — > , < k(x) < 1 and j(x) — > — oo for x — > oo, 

where x is a radial coordinate and a prime denotes differentiation with respect to x. Therefore any 
properties of Eqs. (l)-(3) apply equally well to both the gravitational and electromagnetic geon 
case. Any solutions to (l)-(3) are necessarily equilibrium solutions since the background metric is 
assumed static. In Sec. II. A, the numerical solutions presented in Ref. 1 will be re-examined. The 
results suggest further investigation of the numerical solutions is required in order to determine if 
the boundary conditions are satisfied. In Ref. 6, it was shown that admissible equilibrium solutions 
to Eqs. (l)-(4) exist. However, the stability of these equilibrium solutions were not studied. The 
investigation presented in Sec. II. B constructs a phase portrait of the ODEs from which both the 
analytic form at spatial infinity and the stability with respect to perturbations of the amplitude 
eigenvalues of any solutions will be obtained. Knowledge of solution stability provides a basis for 
investigating the evolution in time of these solutions. Unlike other investigations, 1,7 we apply a 
small amplitude time-dependent perturbation to an equilibrium solution of (l)-(3) for the case 
of an electromagnetic geon. This is done in Sec. III. Solving the time-dependent perturbation 
equations leads to a contradiction with one of the initial assumptions. The contradiction suggests 
that neither an electromagnetic geon nor a gravitational geon is a viable construct since not all of 
the requirements for existence of a geon are met. This interpretation of the results obtained from 
this investigation will be discussed in Sec. IV. 
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A. Numerical integration 



Wheeler 1 originally solved the system (l)-(3) by numerical methods in 1955. Since then computer 
algorithms have evolved considerably for solving differential equations. It is therefore worthwhile to 
utilize modern techniques 8 in re-examining those solutions. Even with the algorithm employed in 
Ref. 1 for solving the equations, Wheeler's results are remarkably accurate. 

The geon problem (both electromagnetic and gravitational) is reduced to finding a solution to 
the autonomous system (l)-(3). Admissible solutions to Eqs. (l)-(3) are defined as those <j>(x), j(x) 
and k{x) that satisfy the following criteria: 

1. For large negative x: The wave function (j)(x) — ► and metric function k(x) — > 1. Under these 
conditions j'(x) — > 2. If <j>(x), j(x) and k(x) are solutions to the autonomous system (l)-(3), 
then so are (f>(x + a), j(x + a) and k(x + a) where a is a constant. Choosing the integration 
constant for j(x) to be zero fixes a and consequently j(x) — > 2x. This removes any ambiguity 
in the start of the integration process. Thus for large negative x, <j>{x) satisfies the equation 



2. For large positive x: It is required that <fi(x) — > 0, < k(x) < 1 and j(x) approach large 
negative values. 

The only free parameter is the amplitude A of the wave and this must be chosen so that the 
solution fits the boundary conditions. The nonlinearity of the problem makes it necessary to integrate 
the system of equations numerically. The integration is started at x = —4. The initial conditions 
are as follows: 




(5) 



The approximate solution to (5) as given in Ref. 1 is 




(6) 



0(-4) - 0o 



(7) 




(8) 



fc(-4) = 1, 



(9) 



i(-4) = -8. 



(10) 
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where <p is to be chosen to give an admissible solution. Those 0o values which yield admissible 
solutions will be referred to as eigenvalues of the system (l)-(3) with initial conditions (7)-(10). 

The behaviour for <j>{x) as x — ► oo depends upon the value of 0o (which translates into an initial 
choice of the amplitude A). Figure 1 illustrates numerically integrated solutions for values of O 
given in table I. Solution set 1 shows that for sufficiently large values of 4>o, <f>{x) reaches a positive 
minimum and then increases exponentially. As O is allowed to decrease, the exponential growth 
of (j)(x) is delayed. This is depicted by solution sets 2 and 3. A further reduction in O results in 
<f>(x) — > — oo exponentially as shown in solution sets 4-6. A possible admissible solution lies between 
solution sets 3 and 4. 

The mass of the geon inside radius p is related to the function k(x) in the following way: 

M(p(x)) = iA (x) = l(l-fc 2 ), (11) 

with b = l/Qi(oo) = fc(oo). This implies that 

< k 2 (x) < 1 as x -> oo (12) 

in order to have a positive total mass. The mass factor k(x) gives a positive mass throughout the 
integrable region and appears to have a fe(oo) value of approximately 1/3. The 'active region' can 
be identified in the x coordinate system as the region where the function j(x) is positive. In this 
region, the function <p(x) has oscillatory behaviour. The function j(x) is positive only for a limited 
range in the neighbourhood of x = 1, thus identifying the active region. In Fig. 1, the active region 
begins at x ~ 0.12 and ends at x ~ 2.13. The first admissible solution (characterized by 4>{x) having 
one local maxima and no local minima) appears to lie between those values of 0o in the range 
9.7904 x 10~ 5 < 0o < 9.7906 x 10~ 5 . Qualitatively, these results are similar to those in Rcf. 1. The 
only main difference between the calculation of Rcf. 1 and the present one is that in Ref. 1, the first 
admissible solution appears to lie in the range 1.03000 x 10~ 4 < 0o < 1.03125 x 10~ 4 and the active 
region starts at x ~ 4.05 and ends at x ~ 6.02. 

B. Existence and stability of equilibrium states 

We are interested in determining the analytical behaviour of the solutions shown in Fig. 1 as x — > oo. 
By constructing the phase portrait for the differential equations, it will be possible to determine 
both the existence and stability properties of potential admissible solutions. We start by rewriting 
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Fig. 1. Results of the numerical integration for the geon differential equations (l)-(3). The values 
of <po for solution sets 1-6 are summarized in table I. The integration started at x = —4 and could 
not proceed beyond x ~ 7 for the given initial values. The active region begins at x ~ 0.12 and 
ends at x ~ 2.13. A possible admissible solution lies between sets 3 and 4. Note that for j(x), 
curves 1-6 are indistinguishable. 



Table I. Values of cf>0 for solution sets 1-6. 

Solution set <po 

1 9.7910 X 10" 5 

2 9.7908 x 10 -5 

3 9.7906 x 10" 5 

4 9.7904 x 10~ 5 

5 9.7902 x 10" 5 

6 9.7900 x 10~ 5 
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equations (1)~(3) as the set of first order equations 

u' = -jk<j>, (13) 

<P' = «, (14) 

k> = -4>\ (15) 

f = 3- k- 2 +u 2 k- 2 . (16) 



It is sufficient to look for solutions with the properties 

<p, u — > 
fc-» 1 

4>, u —* o 

fc — > constant > 



as a; — > — oo, 

(17) 

as x — > +oo 



and j remains finite for finite x. In a phase space, the critical points (or equilibrium points) are 
characterized by those points where the derivatives of u, (f>, k and j are zero. The analytic behaviour 
of the solution about a critical point is determined by analyzing the corresponding linear system in 
a neighbourhood of that critical point. 9 

The first step is to obtain the critical points of the system (13)-(16). One can easily verify that 

u = 0, (18) 
= 0, (19) 
fc = ± -L (20) 

is sufficient to satisfy u 1 = cf)' = k' = j' = 0. Thus the coordinates of the critical point in the phase 
space are 

u = 0, (21) 

= 0, (22) 

fc =-L, (23) 

3 =s, se R, (24) 

where s is any value of j(x). The function k{x) cannot pass through zero since Eq. (3) becomes 
singular. The positive root of Eq. (20) is chosen to ensure positive k(x), since initially fc(— oo) = 1. 
It is useful to shift the critical point to the origin using the following transformation: 

4>{x) = fi{x) i u(x) = f 2 (x), 

1 (25) 
K x ) = h{x) + -j=, j(x) = f 4 (x) + s. 
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Therefore the field equations become 



f[ = 1 2 
= — /l) 



(/ 4 + S )(/ 3+ _J/ 1 , 



/ 4 = 3-(l + /|)(/ 3 + -^ 



(26) 
(27) 
(28) 
(29) 



with the critical point at fx = f 2 = = fi = 0. To linearize the field equations about this critical 
point, a MacLaurin series of = (/i, /2, h, fi) , i = 1, • • • , 4 is taken to first order and evaluated 
at the critical point (denoted c.p. below), i.e. 



fi=f< 



c.p. dfi 



f + dn 



c.p 



h + • ■ 

c.p. 



Written in matrix form, the linearized field equations are 

dw IV /I 

-— = M w, 

ax 

where 



(30) 



(31) 



h 

h 



and M 



/ 1 0\ 

-s/y/3 



V 6V3 Oj 



(32) 



The general solution to the above matrix differential equation is the eigenvector 
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+ C 3 





e aa: + c4 
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V o J 










W 



(33) 



where a 
is 



-s/V3 and a, i = 1, ... ,4 are constants. Therefore the solution to the linear system 



/i(ar)= C3 e Q3; + C4e- Q3: , 
f2{x) — c 3 a e" x ~ C4 a e~ 

h{X) = c 2 , 
/ 4 (a;) = ci + c 2 x. 



(34) 
(35) 

(36) 
(37) 
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The eigenvector w shows that in a neighbourhood of the critical point, the nonlinear system decouples 
into the disjoint subspaces (/i, / 2 ) and (/ 3 , f±). 

The phase space projection of (/ 3 , / 4 ) (which is proportional to (k, j)) is illustrated in Fig. 2. 
Since s can take any value of j, the critical point lies at an arbitrary position s < max(j) on the curve 
k(x) = -j^ c 2 + -j= which is transformed back to the origin in the (/ 3 , fa) subspace. Equations (36) 
and (37) show that in a neighbourhood of the critical point, but not on the /4-axis, / 3 and / 4 behave 
as a constant and a linear function of x respectively. If one is on the / 4 -axis in a neighbourhood 
of the critical point, then c 2 = 0. Hence / 4 = c\ defines a continuous set of critical points. These 
critical points are examples of non-isolated critical points. The (/ 3 , f±) projection is insufficient for 
determining the existence and stability of admissible solutions since it only gives information about 
the functions j and k. 

The nature of the (/1, / 2 ) subspace 10 (or the (<j>, u) subspace) depends upon the parameter s e j. 
It will determine the stability properties with respect to perturbations of the amplitude eigenvalues 
of any admissible solutions. However, a 3-dimensional phase space projection in the coordinates 
(<j), u, j) is necessary to determine the existence of admissible solutions. An illustration of the two 
possible phase space projections of (</>, u) are shown in Fig. 3. Examining (34) and (35), if s > 0, 
then <f> and u behave as sinusoidal functions of x. This type of critical point is described as a center 
(Fig. 3(a)) and is a stable critical point. If s < 0, then <f> and u have an exponential behaviour and 
the critical point is unstable. This type of critical point is described as a saddle point (Fig. 3(b)). 
By following the integration procedure in the parameter x for solution sets 1 and 6 of Fig. 1, the 
behaviour of the complete nonlinear system can be described. 

Figure 4 shows the 3-dimensional ((j), u, j) phase space projection of solution sets 1 and 6. The 
integration procedure starts in plane B of Fig. 4 at j = —8. In addition, the solution trajectories 
start somewhere along the line u = (1/16 + v / 8)0o which must necessarily lie to the left of the 
unstable asymptote £+. One such point is labelled "1" in Fig. 3(b). In order for the system to 
satisfy the boundary conditions (17), it is necessary for at least one solution to flow along the 
unstable asymptote £~ (in the j = s — * —00 plane). If Fig. 3(b) were a complete description of the 
phase space, then it would be impossible for a solution starting at position "1" to cross £ + . This 
is a consequence of the well-known property of autonomous systems that phase space trajectories 
do not cross. However, as the integration process in x continues, the value of j increases from a 
negative value to a positive value. Therefore the nature of the critical point in the 2-dimensional 
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Fig. 2. The (/3, fn) phase space projection illustrates a non-isolated critical point. The critical 
point under investigation is located at the origin. The /4-axis is a continuous set of critical points 
in a neighbourhood of the origin. The functions /4 oc j and fsock behave as a linear function of 
x and a constant function respectively in a neighbourhood of the critical point not on the /4-axis. 
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Fig. 3. (a) Illustration of a stable critical point. This type of critical point occurs when the value 
of the parameter s > 0. (b) Illustration of an unstable critical point. This type of critical point 
occurs when the value of the parameter s < 0. 
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Fig. 4. The 3-dimensional (<fr, u, j) projection of the phase space curves for the numerical solution 
sets 1 and 6 shown in Fig. 1. The numerical integration starts in plane B, proceeds up through 
plane A and continues down through plane C. The nature of the critical points along the j-axis 
(origin of the (</>, u) planes) change from centers (j > 0) to saddle points (j < 0) demonstrating 
the system instability. 
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((f), u) phase planes changes temporarily from a saddle point to that of a center. Plane A of Fig. 4 
is an illustration of one such critical point. Soon afterward, j decreases to negative values and the 
critical points are saddle points once again. However, the solution trajectories have crossed £ + and 
now follow the flow along the unstable asymptote £~. Upon the transition of the critical point 
from centers to saddle points, the asymptotes have been reestablished with solution sets 1 and 6 on 
opposite sides of the £~ asymptote. The positions where solution sets 1 and 6 cut the (4>, u) planes 
for j < are schematically illustrated as points "2" and "3" respectively in Fig. 3(b). In Fig. 4, 
these positions are most clearly seen in plane C. Since trajectories for autonomous systems do not 
cross and the trajectories depend continuously on the initial data, there must be a value of <fi for 
which the trajectory approaches £~ as j = s — > — oo. The existence of this trajectory shows that it 
is possible to find an eigenvalue of 4> which satisfies the boundary conditions (17). However, the 
nature of the critical point as j = s — ► — oo (x — ► oo) is a saddle point and therefore this eigenvalue 
solution is an unstable solution. Since solution set 1 cuts plane C at position "2" of Fig. 3(b), 
the flow of the integration process requires this set to approach . Similarly, solution set 6 must 
approach £~ , since it cuts plane C at position "3." Hence, any small perturbation of the eigenvalue 
of <p implies the constant c 3 ^ in (34) and (35). Thus the non-eigenvalue solutions do not satisfy 
the boundary conditions. Figure 5 shows the (</>, u) subspace for the six solution sets of Fig. 1. The 
projection is for j < — 15. A comparison of Fig. 5 to Fig. 3(b) confirms that the admissible solution 
is unstable. 

The existence of admissible solutions and instability of the electromagnetic geon system were 
discussed in Ref. 1. However, it was based on the numerical solution curves similar to Fig. 1. 
Performing a phase portrait analysis, we have formally shown that an eigenvalue does exist 11 which 
satisfies conditions (17). We have also shown that this solution must necessarily be unstable with 
respect to perturbations of the eigenvalue of fa. This result suggests that geon constructs are 
dynamically unstable, i.e. the ensemble of waves must collapse or explode. In Ref. 1, it was also 
suggested that a spherical geon would most likely collapse to a toroidal geon, 12 presumably thought 
to be more stable. It was argued by Ernst 12 that the construction of a toroidal geon could be 
realized if one had complete knowledge of a linear geon (which approximates a small segment of the 
toroidal geon). Numerical evidence for amplitude eigenvalues in analogy with the spherical geon were 
presented in Ref. 12. This work was extended in Ref. 13 by performing a phase portrait analysis. 
It was found that only unstable admissible solutions exist, as in the case of the spherical geon. 
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Fig. 5. The (0, u) phase space projection of solution sets 1-6 for j < —15. It has the characteristics 
of Fig. 3(b). This indicates that the critical point at the origin is unstable as j — > — oo. 



15 



It has been suggested 6 that this is sufficient for proclaiming the existence of both types of geons 
(electromagnetic and gravitational). However, in Sec. I, it was noted that a true geon must have the 
property that the ensemble of waves comprising the geon be confined on a time-scale much longer 
than the typical period of the constituent waves. Otherwise it would not be possible to attribute a 
structural form to the geon. What is yet to be determined analytically is the dynamic behaviour 
of geon constructs when perturbations are present. The next section presents a time-dependent 
perturbation analysis of the equations that describe the electromagnetic geon in an attempt to 
determine this time-scale. 

III. Time-Evolution Analysis of the Electromagnetic Geon 

The electromagnetic geon model employed in Ref. 1 assumed a background metric independent of 
time. This precludes the possibility of studying the time-evolution of the individual modes coupled 
to a time-evolving background metric. Instead of introducing a time-dependence into the system 
of equations and solving the coupled wave-background system, an alternate method was employed 
in Ref. 1 for determining the time-scale of collapse or 'lifetime' of the electromagnetic geon. The 
approach was to study the behaviour of a single photon (borrowing particle physics terminology) in 
an effective potential created by the ensemble of photons comprising the geon. It was expected that 
photons would leak out of the potential well in an irreversible dissipation of energy. The model for 
the rate of photon leakage was based on the quantum mechanical process of radioactive alpha-decay. 
From this model, the probability of barrier penetration, a, (called the 'attrition' in Ref. 1) has a 
value of the order 

a~e- 2P , (38) 

where P is the barrier penetration integral (c.f. the Gammow factor of alpha-decay 14 ). The function 
P is roughly proportional to the height of the barrier as encountered by the escaping particle. A 
small attrition would translate to a small rate of radiation leakage out of the effective potential. The 
time to collapse or the lifetime of the geon in this model is inversely proportional to the product of 
the attrition and angular frequency of the vibrational mode under consideration. It was estimated 



using data from the numerical integration performed in Ref. 1 that the attrition a ~ e _1 ' 52 V / ^+ 1 ^. 
As the angular momentum I is increased, the attrition decreases exponentially which suppresses the 
photon leakage. Thus it was concluded that it was possible to construct a geon of virtually infinite 
lifetime. 
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The only similarity between the concept of a geon leaking photons and radioactive alpha-decay 
is the shape of the potential well found in each problem. However, there are important differences 
to consider. It is to be stressed that a geon is a classical object built from classical (massless) 
fields. Radioactive alpha-decay is a Coulomb repulsion effect involving massive, electrically charged 
particles. The mechanism for the emission of an alpha particle from a nucleus is the quantum 
mechanical effect of 'tunnelling' through a potential barrier. In a classical system, it is not possible 
for a particle to tunnel though a potential barrier. In addition, the concept of photon leakage 
implies the background metric (or potential barrier) evolves on some time-scale. An alpha-decay 
model cannot describe such an evolution. The quantum mechanical nature of alpha-decay brings into 
question the validity of using such a model for determining the lifetime of a classical object such as a 
geon. There are known phenomena which represent classical wave penetration of a potential barrier. 
For example, the optical phenomenon of frustrated total internal reflection 15 is such a process. It 
is important to note that in both alpha-decay and frustrated total internal reflection, the potential 
is supplied by some material and not the waves themselves, as is necessary for the case of geons. 
Whether the analogy exists between these examples of barrier penetration and the time-evolution 
of a geon based upon the coupled Einstein-Maxwell equations is the subject of this section. 

Another approach towards determining the lifetime of an electromagnetic geon is that of Brill. 7 
The method is to study the evolution of the ensemble of photons which produce the effective potential 
using a thin-shell model for the electromagnetic geon. It was found that the radial position of the 
thin shell underwent a displacement towards collapse. It was also stated that the rate of collapse 
was 'slow.' 

The junction condition problems associated with analyzing thin-shell geon models has previously 
been discussed in some detail in Ref. 5. In addition, the evolution of the thin-shell model in Ref. 7 
does not allow for leakage of radiation during the collapse nor does it allow an evolving shell thickness 
(evolving active region) as one might expect. The effect of correcting these deficiencies on the rate 
of collapse is not clear. A full understanding of the evolution of a geon must take into account 
the evolution of the typical individual modes of vibration coupled to the evolution of the collective 
ensemble of waves in a singularity-free model. 

It is evident that models for studying the evolution of geons must be based on solving the 
Einstein (or Einstein-Maxwell) field equations. Intuitive models arc not sufficient for describing the 
true physical system. To avoid the interpretation problems associated with the alpha-decay and thin- 
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shell models and correct for the deficiencies of each model, the derivation of the electromagnetic geon 
equations 1 will be modified to permit the study of the time-evolution of the electromagnetic geon. 
The electromagnetic geon will be studied instead of the gravitational geon because of the relative ease 
in computation for the former. Equilibrium solutions for the gravitational and electromagnetic geon 
are governed by the same set of ODE's. Therefore it is not expected that the modified gravitational 
geon equations in the high-frequency approximation would yield a significantly different result from 
the electromagnetic case. 

We are interested in following the evolution of the electromagnetic geon in the radial direction. 
Observing the time-evolution of the metric reflects the evolution of the ensemble of electromagnetic 
waves comprising the geon. However, it is not sufficient to simply perturb the background metric 
functions, ft is the electromagnetic waves which are the source for the gravitational field, hence both 
the gravitational and electromagnetic quantities must be perturbed. This will be done by applying 
an amplitude perturbation on the electromagnetic waves comprising the geon in such a way as 
to induce the background metric to evolve in time. Frequency perturbations are not explicitly 
considered in the following derivation for two main reasons. Firstly, the stability analysis of the 
previous section indicates that the instability of the admissible equilibrium solution originates from 
changes in the amplitude eigenvalue (initial condition) . Secondly, it can be shown that a perturbation 
of the form ft — > Q + SCI is a special case of the amplitude perturbation studied below. Restricting 
study to the radial direction maintains the field equations in their simplest form. It should be 
emphasized that the time-space average of the electromagnetic disturbance must be incorporated 
into the background metric equations to maintain spherical symmetry but still allow for the solution 
to evolve in time. This time-averaging problem will be addressed when the perturbation is applied. 
Before the perturbation analysis is performed, the angle-averaged time-dependent electromagnetic 
geon field equations must first be derived. This part of the derivation follows closely that of Ref. 1. 

The equations presented below are derived in greater detail in appendix A. Only an outline of 
the derivation is given here. We start by defining the electromagnetic vector potential for one mode 
of the electromagnetic waves 



= (0, 0, 0, A v ) , 



(39) 



where 



A v = a(r, t)B{6) , B(0) = sin 8 ^-Pt (cos 6). 



(40) 
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The function a(r,t) has been left unspecified at this stage. The time-dependent background metric 
is 

ds 2 = g aP dx a dx fi = -e^dt 2 + e A dr 2 + r 2 d9 2 + r 2 sm 2 9d(p 2 , (41) 

where 

v = v(r, t), A = A(r, t). 

In the absence of charges and currents, Maxwell's equations in a curved space-time are 

1 d 



{sf—gF 1301 ) = 0, (42) 



\f— ~g dx° 

F a p n + F iaJj + F M . a = 0, (43) 

where g is the determinant of the metric (41) and the Maxwell tensor, F a/ 3 is related to the four- 
vector potential as F a p — Ap ia — A ai p. The Einstein equations for the electromagnetic geon are 

G/=8tt(T/), (44) 

where ( • ) denotes a time-space average over all N active modes 16 of the electromagnetic waves. 
Substituting (39) into (42) and (44), taking the angle average and finally transforming to the p = Qr 
coordinate system yields the wave equation 

o d 2 a o*o of 2L\ o ( 2L\ 2 d 2 a , s 

n V - 1 p vt) q vt) Q d^ =0 (45) 



and the background field equations 



+ fi 2 rV 2 (l-y)Q(a 2 ) T ^ (46) 



dQ nf ( 2 / ( 9a \ 2 \ / ( da N 



dp* p—2L\ \\dp* ) J \\dp 



9L n 2 n 2 Q^/^\ 



^ ^ + ^ . (47) 



at-'-^WW; (48) 
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where 



Ot* 2 

and 



(-^"'^(O-t)"^!)' 



Nl* 2 



K ^V27TT r = V*(* + i)- (si) 

In the above equations L{p) and Q{p) are metric functions (see Eqs. (A24), (A25)) and the symbol 
( • ) T denotes a time-average is to be taken. Equations (45)-(48) are the starting point for devel- 
oping the dynamic perturbation (time-evolution) equations. The dQ/dt equation (found from the 
G g e — 8n (T g 9 ) equation) is not used in the subsequent analysis, but is given in appendix A for 
completeness. 

The time-averaged equilibrium solution in Ref. 1 has the form 

ki f2 a(p, t) = fo(p) sin O t, (52) 
L(p, t) = L (p), (53) 
Q(p, t) = Qo(p). (54) 

where fo(p), L (p) and Qo(p) are known functions. 17 We will designate this as the 'unperturbed so- 
lution.' Two general forms for the radial perturbation of the wave function a(p, t) will be considered. 
The first is given by 

Kl n a(p, t) = , f {p) sinn t + Su 1 (p,t) + S 2 u 2 (p,t) + 0(S 3 ), S < 1. (55) 

where 5 is the expansion parameter. (Note that the addition of a phase constant to sin Sit does 
not affect the results which follow. Thus the phase constant is set to zero.) As a result, a small 
time-dependent perturbation is introduced in the metric functions 

L(p, t) - L (p) + SL^p, t) + S 2 L 2 (p, t) + (S 3 ) , (56) 
Q(P, t) - Qo(p) + SQ^p, t) + S 2 Q 2 (p, t) + (S 3 ) . (57) 

The perturbation expansion will be carried out to the first order in S. Before the perturbed system 
is solved in a self-consistent manner, the problem of time-averaging the functions on the right hand 
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side of Eqs. (46)~(48) must be addressed. From the definition of a time-average and Eq. (55), 



where T = 27r£l~ 1 is the period of the electromagnetic waves. In order to develop the perturbation 
analysis, it is necessary to make some assumptions about the function ui(p,t). The presence of 
unevaluated integrals on the right hand side of the differential equations (46)-(48) would not make 
it possible to proceed with the analysis beyond this point. Let us suppose the time-dependence of 
u\(p,t) was sinusoidal and its characteristic frequency was of the order £1 of the electromagnetic 
waves. In this case, the time-dependence is lost to all orders in <5 upon time-averaging. In essence, 
this assumption on u\(p, t) precludes the possibility of a time-dependent evolution of the system. 
This is not satisfactory. To maintain a time-dependence after time-averaging, another time-scale 
will be introduced into the problem. Suppose the time-dependence of Ui(p,t) was again sinusoidal 
and its characteristic frequency was of the order oj -C O. In this case, u\(p,t) is approximately 
constant over the short time period T = 2ttQ~ 1 of the electromagnetic waves and thus is constant 
in the time-average integral. Evaluating the time-average of a 2 (p, t) under this assumption yields 



Therefore the time-dependence is not present until the second order in 5. This is sufficient to proceed 
with the time-evolution of the system, since each order in the expansion parameter 5 must be set to 





-f^p) + 6 2 ul(p, t) + 0(S 3 ). 



(58) 
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zero. To first order in S 



K 2 ft 2 (a 2 (p,i)> T = i/ 2 + O(5 2 ). (59) 



2 

Similarly, the remaining time-averages are 

»?f 2 ((|) 2 ) = ^V„ 2 + O , (61) 

Substitution of (55)-(62) into (45)-(48), expanding to first order in 5 and setting each order in 6 to 
zero yields the unperturbed equations (63)-(65). The properties of the unperturbed equations 



dp 



~ 'i_r 2 Q 2 p- 2 (i- — ) )./„ = (). m) 



d fo , I i 7*2^,2 -2 A 2i 



^0 1 ( j-2 , ( d/o\ 2 , ,*2^,2 -2 f-, 2L \ • 



are known from Rcf. 1 and therefore can be used in the analysis of the first order equations. Setting 
the first order part of the wave equation (45) to zero yields 

A(p, t) sin n t + B(p, t) cos fit + C(p, t) = 0, (66) 



where 



A(p, t) = Sll (Qi(p, t) (l - - 2p~ 1 L 1 (p, t)Q ) x 



dp V dp ) dp 2 V 9 J dp dp V P 
+ 2p : 7 -V) (l /) ? ;0 / ) ( ] _^£^/ + Q (l-^ 



x ^ 2p ~ d^ (v Ql (v L ° " p ~^) + Qo [ LliPl t] ~ p dp~ Ll{p ' t] 

(( 2L \ d , dQ \ d/ 
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+ (Qi(p, t) (l - ^) - 2p- 1 ii(p, t)Qo) 0j 



■ 1 ~^y) Qo 1 (Qi(p,t)(^-^)-2 n - > Ld/>.l)Qn)j\>). (07) 



and 



c( P , t) ee ^j«i(p, t) - r 2 Q§p- 2 ^1 - Ul ( P , t) n- 2 — Ul ( P , t). 



(68) 



(69) 



In order to satisfy the first order wave equation (66), the conditions A(p, t) = 0, B(p, t) = and 
C{p, t) = must be imposed. This is justified since sinf2t and cosfH are independent functions in 
the approximation where A(p, t), B(p, t) and C(p, t) are slowly varying functions of time. We will 
focus upon the latter equation, since it is the simplest of the three equations. Setting (69) to zero 
yields the equation 

dp* 

To solve this equation, we will first use the method of separation of variables. Later, this restriction 



|2 / 2L \ d 2 

^ Ul ( P , t) - r 2 Ql P - 2 ( i - -^J ui{p, t) - n- 2 — Ul ( P , t) = o. 



(70) 



will be removed. Let 

Ul (p, t) = u{p)T{t). (71) 
Substituting (71) into (70) and dividing by u{p)T{t) yields the two ordinary differential equations 

d2?i(p) - 1* 2 Q 2 P- 2 (l- 2 -^) u(p) + 0u(p) = 0, (72) 



dp* 2 V p 



-^±l +/ 3n 2 T(t) = 0, (73) 

where —0 is the separation constant. The solution to (73) is 

T(t) — c 3 sin(u;t + c 4 ) , C3, C4 - constants, (74) 

where we have chosen = uj 2 /Q 2 , Making this choice for satisfies the requirement of (58) 

for a meaningful time-average. 
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In order to solve Eq. (72), it is necessary to apply the high angular momentum approximation. It 
is therefore necessary to transform C (p, t) to the x coordinate system and expand in inverse powers 
of P 1 / 3 as is done for the unperturbed system. 1 The transformation in Ref. 1 for the unperturbed 
functions is repeated here for convenience 

x = (p* -P)?*- 1 / 3 , (75) 
dp* = l* 1/3 dx, 

P = r + r 1 /\ (x) + -- - , 

Lo = l*Xo(x) + l* 2/3 Xi(x) + l* 1 / 3 X 2 (x) + • • • , (76) 
Qo = l/k(x) + l*- 1/3 qi (x) + l*- 2/3 q 2 (x) + ■■■ , 

fo = r 1 ' 3 ^) + 4>i (x) + r-v z <h{x) + ■■■ . 

In addition, the function u(p) must be expanded in a similar manner, i.e. 

u = l* 1/3 p (x) + p^x) + l*~ 1/3 P2{x) + ■■■ ■ (77) 
After a lengthy computation, the asymptotic expansion of (72) yields 

In order for (78) to be satisfied for large arbitrary I* (in the limit I* — > oo), each order of P 1 / 3 must 
be set to zero. Since setting po( x ) = implies the absence of electromagnetic wave perturbations, 
the bracketed term must be zero. It is known from the unperturbed system that 1 

\ ) {x) = \(l-k 2 {x)). (79) 

Substituting (79) into (78) and setting the bracketed term to zero yields the relation 

lu 2 = n 2 (so) 

which must hold in order to satisfy the condition C(p, t) = 0. However, Eq. (80) is a contradiction 
to the original assumption ui <C fi. Because of the presence of the contradiction, it is not necessary 
to solve the remaining field equations. 

We arrived at this contradiction through the assumption that (70) could be solved by separating 
variables (Eq. (71)). The same result is obtained if u\(p, t) is not separable as is shown below. 
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Modifying (77) as follows 

«i = l* 1/3 Vo(x, t) + m{x, t) + l*- 1 / 3 ^, t) + • • • . (81) 
and expanding (70) in inverse powers of I* 1 / 3 yields 

-rV3 (n-'gy m*> *) - 1 p 2 ( t°) (x) ) t) + o(i) = o (82) 

to lowest order in Z* 1 / 3 . Setting each order in I* 1 / 3 to zero requires {j,q(x, t) to satisfy the differential 
equation (using (79)) 

d 2 

-^ IM) {x,t)+n 2 vo(x,t) = 0. (83) 

The solution is 

fj, (x, t) = c 5 (x) sm(Qt + c 6 (x)). (84) 

The characteristic frequency of ^0(^7 t) is f2 which contradicts the assumption u\(p, t) <~ /Uo(a;, i) ~ 
uo < 0. 

Before we proceed with the discussion of the perturbation analysis based on (55)-(57), a second 
form for the wave function perturbation should be considered. Consider the perturbation in the 
form 

ni Qa(p, t) = (fo(p) + 6 Ul (p, t)) sin fit + O (<5 2 ) , (85) 

where the characteristic frequency of Ui(p, t) is of order u <C ft. This form can be interpreted as 
a slowly evolving amplitude of the rapidly varying function sinQi. Equations (55) and (85), in 
addition to the assumptions placed on ui(p, t), cover the entire range of possibilities for these types 
of perturbations (for example, a perturbation of the form a(p, t) — > a(p + S£(p, t), t) reduces to 



(85)). Evaluation of the time- averages yields 

k\ n 2 (a 2 ( P , t)) T = \f 2 + 6f oUl + O (S 2 ) , (86) 

"'""((S) 2 ) =j tf /„ 2 + «n 2 /o«i + o(« 2 ), (88) 
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Hence, the time-dependence of the time-averaged functions becomes manifest at the first order in 
5. This greatly increases the mathematical complexity of the system. The detailed analysis for 
this system is derived in appendix B for the case of Ui(p, t) separable. From this analysis, the same 
contradiction results for the characteristic frequency of Ui(p, t) as that found from the analysis based 
upon Eq. (55). For nonscparablc u\{p, t), the differential equations become unmanageable. Hence, 
it was not possible to obtain a conclusive result. However, due to the similar nature of the systems 
based upon Eqs. (55) and (85), there is no reason to suspect a different result for the nonseparablc 
case. 

IV. Discussion 

We have seen from the analysis of the static solution (Sec. II) that the confinement of gcon constructs 
demands an absolutely critical choice of initial condition (amplitude eigenvalue). The slightest de- 
viation from that choice leads to a totally unconfincd structure. While some might argue that the 
confined structure being indicated with the critical initial condition is already satisfactory, 6 the gen- 
eral experience with solitonic structures is one of essentially confined solutions in the neighbourhood 
of the best choice of critical condition. 18 The failure to find a family of near-confinement in the case 
of the geon already raises suspicions as to its viability. 

In the previous section, the time-evolution of an electromagnetic geon equilibrium solution was 
studied with the objective of determining the time-scale of collapse away from the equilibrium 
configuration. Perturbations of the form (55) and (85) were analyzed under certain assumptions. 
The problem of time-averaging the source terms (right hand side) of the field equations (46)— (48) 
requires the characteristic frequency u> of the perturbation term of the wave function u\(p, t) to be 
much less than the characteristic frequency of the unperturbed solution. Without this assumption, 
the time-dependence of the perturbations is lost upon time-averaging, to all orders in the expansion. 
Hence, the assumption u\(p, t) ~ lu ~ Q is not satisfactory for studying the time-evolution of gcons. 

The assumption u\(p, t) ~ u <C Q, for both forms of the perturbation, solves the time-averaging 
problem in a simple manner. With this assumption, the differential equations maintain a time- 
dependence after the time-average has been taken over the typical period of the high-frequency 
waves. The perturbation analysis leads to the requirement u> ~ Q in order for the field equations 
to be satisfied. This is a contradiction to the original assumption u -C O. Since all the possible 
combinations for the form of the perturbation (and assumptions placed on U\{p, t)) have been 
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explored, the possible interpretations of the results of Sec. Ill are given below. 

A reasonable interpretation of the contradiction in the time-evolution analysis is that the condi- 
tion of slow evolution of the background cannot be satisfied. Since not all of the required conditions 
are satisfied, it is not possible to construct a geon comprised of high-frequency waves. It could be 
argued that an electromagnetic geon could be built from low-frequency waves. This has not been 
ruled out by our model since it only accommodates high-frequency waves. However, a gravitational 
geon necessarily must be constructed from high-frequency waves, for otherwise the effective stress- 
energy would not be of the correct order of magnitude to create the background gravitational field 
binding the waves. 3-5 Therefore the gravitational geon studied in this paper is subject to the same 
fate as its high-frequency electromagnetic counterpart. 

A geon with a rapidly evolving background metric, where the background is somehow regarded 
as being distinct from the small amplitude waves on the background, is conceptually unsound. The 
definition presented in Sec. I for this type of geon requires the background solution of the Einstein 
or Einstein-Maxwell field equations be quasi-stable on a time-scale much longer than the period of 
the constituent waves. If the background metric evolves away from the equilibrium configuration on 
the same time-scale as the constituent waves, one cannot speak of the waves binding gravitationally. 
Under these circumstances, there is no geon structure to identify. Apart from this argument, the 
realization of a geon with a rapidly varying background metric 7^ is problematic for another reason. 
As discussed in Ref. 5, if a spherically symmetric background is allowed to vary harmonically with 
frequency Q, comparable to the frequency of the gravitational or electromagnetic waves, one expects 
a parametric resonance 19 for the modes with u n = nCl/2, with neN. The strength of the resonance 
is a maximum for n = 1 and decreases rapidly as n increases. In the limit of a static background, 
the resonance phenomenon disappears. Accordingly, on the basis of studies of perturbations of black 
holes and relativistic stars, 20 it is expected that in the case of a stationary axisymmetric background 
metric describing a rapidly rotating geon, the resonance phenomenon between the perturbations and 
the background metric occurs. In the general case of a time-dependent and rapidly varying back- 
ground metric 7^ (t, x) without symmetries, it is not known how to decompose metric perturbations 
on a complete set playing the role of the tensor spherical harmonics in the spherical case, or even 
how to define frequencies in the strong curvature region. However, if such concepts can be given a 
meaning, it seems reasonable to expect some kind of resonance phenomena between the background 
metric and its wave perturbations. All these resonance phenomena certainly do not contribute to 
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the realization of a stable configuration, but rather are associated with instabilities that tend to 
disrupt the system. 

The contradiction which arises in the perturbation analysis of Sec. Ill is interpreted as a break- 
down of the model, i.e. the perturbation model is not able to describe the evolution of the physical 
system. It might be argued that an alternative method of implementing the time-evolution (not using 
perturbative methods) may be better suited to determine the time-scale of evolution. For example, 
it may be possible to develop an exact numerical solution of the full Einstein (or Einstein-Maxwell) 
equations without the splitting of the metric into a background and waves on the background or tak- 
ing time-averages. Considering the complexity involved in analyzing the simple perturbative model 
of Sec. Ill, any new model would undoubtedly be more complex to solve. However, it should be pos- 
sible to approximate any exact method with an appropriate perturbation expansion (as is presented 
in Sec. III). Therefore, it is appropriate to consider the contradiction in a physical context, as was 
discussed earlier. 

At this point, we recall the original motivation which led the authors to re-open the issue of geons 
and their viability. One of the authors 22-24 had been led to propose a new hypothesis regarding 
the localization of energy in general relativity, namely that energy was localized in regions of non- 
vanishing energy- momentum tensor T^ v . There were various factors leading to this. Firstly, the 
energy-momentum conservation laws 

T» v . v = (90) 

are devoid of content in vacuum, producing the empty identity = 0. However, when (90) is re- 
expressed as a vanishing ordinary divergence to create a global form of the conservation law involving 
the introduction of pseudotensors, it is used to compute supposed fluxes of gravitational field energy 
in vacuum where the originating law is devoid of content. It was proposed that the ambiguity of the 
pseudotensorial flux vectors actually reflects the illegitimate injection of supposed physical content 
where none actually exists. Furthermore, it was shown that for Kerr-Schild metrics, all components 
of the gravitational pseudotensors vanish 25 and gravitational plane waves can be expressed in Kerr- 
Schild form. Since a wave is plane in a relatively small region, this is further support to the belief 
that waves of gravity are not actually carriers of energy in vacuum, in accord with the localization 
hypothesis. Other aspects to support the hypothesis had been outlined including the relationship to 
the important earlier papers of Nissani and Leibowitz, 26 the basis for non-excitation of a Feynman 
detector and the work of Virbhadra 27 which showed that localization of energy is confined to the 
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T^ v regions for static and stationary spacetimes. The gravitational geon remained an outstanding- 
challenge to the localization hypothesis since a purely gravitational non-singular construct displaying 
an unambiguous mass would be a clear counter-example. The present work adds new support for 
the hypothesis apart from the value of understanding this basic construct. From another viewpoint, 
the results which we found in this paper are not surprising. From studies extending over 65 years, 
it was recognized that non-singular soliton structures to model elementary particles are not easily 
achieved and they are successful only with a careful mixture of different types of fields (see Ref. 18 
for a review with earlier references contained therein). The electromagnetic geon depends only on 
the electromagnetic field and its own gravity while the gravitational geon is even more restrictive, 
being a purely gravitational construct. In the light of earlier studies, it is not surprising that such 
simple ingredients should resist compactification. 

V. Conclusions 

The construction of a satisfactory gravitational geon model requires an asymptotically flat, self- 
consistent solution of the Einstein field equations which meets the regularity conditions for a 
singularity-free space-time. Furthermore, it must be demonstrated that the evolution in time of 
the geon must take place on a time-scale much longer than the characteristic period of the con- 
stituent waves (quasi-stability property). 

To satisfy these conditions, it was proposed 3-5 that a satisfactory gravitational geon model must 
be constructed in a manner similar to that of Wheeler's 1 electromagnetic geon. This type of model for 
the gravitational geon is in contrast to the thin-shell model of Brill and Hartlc. 2 In order to construct 
a gravitational geon in principle, it was previously established 3-5 ' 21 that gravitational waves of 
high-frequency were necessary. The application of the high-frequency approximation reduced the 
gravitational and electromagnetic geon problem to the same set of ordinary differential equations 
and boundary conditions. Since the background metric is initially assumed static, any solutions 
are necessarily equilibrium solutions. From a phase portrait analysis of the ordinary differential 
equations governing gravitational and electromagnetic geons, it was possible to determine both 
the existence and stability properties of equilibrium solutions with respect to purturbations of the 
amplitude eigenvalues. It was found that admissible equilibrium solutions were unstable to changes 
in the amplitude eigenvalues. Since a basic requirement for the existence of both types geon is 
the quasi-stability property, an investigation of the time-evolution of an electromagnetic geon was 
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performed. In contrast to other investigations, a small amplitude time-dependent perturbation 
to an equilibrium solution was applied. The time-averaging problem is solved by assuming the 
characteristic frequency of the perturbations vary on a time-scale much longer than that of the 
waves comprising the electromagnetic geon. This is in accordance with the requirement that the 
background metric be quasi-stable. Solving the time-dependent perturbation equations leads to the 
characteristic frequency of the perturbations being of the same order in magnitude as the waves 
comprising the electromagnetic geon. This is a contradiction to the original assumption. Thus it 
could not be shown that the time-evolution of the electromagnetic geon proceeds on a slow time- 
scale using standard perturbation theory modified for time-averaged fields. With not all of the 
requirements for the existence of an electromagnetic geon being satisfied, it cannot be concluded 
that an electromagnetic or a gravitational geon is a viable entity. 

Given the results as applied to the gravitational geon, such a construct cannot be considered a 
counter-example to the energy localization hypothesis as discussed in Ref. 22-24. 
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Appendix A. Time-Dependent Electromagnetic Geon Equations 

We start by defining the electromagnetic vector potential for one mode of the electromagnetic waves 

Afi = (0, 0, 0, A v ) , (Al) 

where 

A v = a(r,t)Bi(9) , B t (6) = sin 9-^ Pi (cos 9). (A2) 
The time-dependent metric is 

ds 2 = g al3 dx a dx fi = -c v dt 2 + c x dr 2 + r 2 d9 2 + r 2 sm 2 9d(p 2 , (A3) 

where 

v = v(r, i), A = A(r, t). 
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In the absence of charges and currents, Maxwell's equations in a curved space-time are 

1 d 



(V~gF pa ) - 0, (A4) 



\f— ~g dx° 

F a/ 3,-f + F ia ,p + Fpy.a = 0, (A5) 

where g is the determinant of the metric (A3) and the Maxwell tensor, F a p is related to the four- 
vector potential as F a p — Ap tCt — A a .p. The only nontrivial equation is for a = <p in (A4). It yields 
the wave equation 

° 2a ^ + 1 U-e^|^=0, (A6) 



dr* 2 r 2 dt* 2 



where 



dr* dr' dr* 2 dr \ dr J 

d_ _ < v -\)/ 2 d_ ^_ _ (y-\)iid_ (jv-\)/2?_\ (M) 

dt*~ dt' dt* 2 m\ dtj- { ' 

The Einstein equations for the electromagnetic geon are 

<V=87r(T/>, (A9) 

where ( • ) denotes a time-space average over all N active modes of the electromagnetic waves. In 
the equations below, the energy-momentum tensor for a single mode of electromagnetic radiation is 
given by 

T w ; y = i ( - \r<#F**i) . ( Ai °) 

with F a p defined above. We will only evaluate the angle average of T " . The time-average will be 
dealt with in the main text (Sec. III). In addition to the three angle averages 1,28 



TA 







T r r ) =^ f X W) sinOdO, (A12) 



2 Jo \ W 

T e ^ = (V) TA = ^£(T (1) +T (i) /) T Sin,d., (A13) 
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there is an additional average which represents the radial flow of energy 



T 



I TA 2 



J (T(i)r*) sin<9d<9. (A14) 



Evaluating (A11)-(A14) using (Al) and the integrals of appendix C one obtains 

r/) .3^^)+^ "W), (Ai 9) 



TA 



8nr 2 (2l + 1) 



Je / TA ^8^(2Z + l)\ a / T ' (M7) 

r -*) Nl lV ] J "ra t ) ■ (A18) 

r /ta Airc v r 2 (2l + 1) \ ' ' 7t v ; 

The components of the left hand side of (A9) are 

'A, r (A19) 
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(A20) 



1 



,tt ~ ^ Af t J J (A21) 

(A22) 

The final step in obtaining the time-dependent electromagnetic geon field equations is to make the 
transformation to the p coordinate system. In addition to the transformation 

r = £, (A23) 

we introduce the two metric functions L{p, t) and Q(p, t) through the defining equations 

1 2L(p, t) 
P 

2 A +" = Q 2 (p, t), (A25) 
Q 2 (p, t). (A26) 



e -A = j / A24 n 

P 



l 2L{p, t) 
P 

The operator has the following form in the p coordinate system 



^■^-"('-tX- ,A27 > 
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By defining 

the operators of (A7) simply transform as 

^-=«^-. #T=^#2- (A29) 
dr* <9p* ' dr* 2 dp* 2 ' 

The operator of (A8) transforms as 

2 f 2LY 1 ,d ( f 2LY 1 ,d\ /A N 



dt* 2 \~ p J "* dt\V p ) "° dt 

After applying the transformation (A23) and Eqs. (A24)-(A30), a lengthy but straightforward com- 
putation yields the wave equation 



and the background field equations 



+ fi 2 rV 2 (l-y)o(« 2 ) T ^ (A32) 



9L 9^9^ i / da da\ , ,„.n 

— = njtfQ- 1 ( — — ) (A34 
dt 1 ^ \dp*dt/ T y ' 



and 



d 2 L , /5L\ 2 ^_ x dLdQ 1/ 2L\ 2 2 



dt 2 1 \dt J ^ dt dt 2 V P 

x (>t(p) + B(p) - 2 K 2 /* 2 ^V _4 (a 2 ) ) • (A35) 

In the above equations we have defined 



^ r ^ v^TTi), (A36) 
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*'^^(«-'f('-7) +2 '-'h-| 



(A37) 



and 



dp 2 \ dp 



+ 2P~' 1 



p ) \dp 
2L\ _1 fdL 



+ f 2 0-^-2,-fl-^V 1 ('^-,- i || I. (A38) 



dp \ p J \dp 

Equations (A31)-(A34) arc Eqs. (45)-(48) of Sec. III. 

Appendix B. Perturbation Analysis of a Slowly Varying Amplitude 

To investigate the time-evolution of the unperturbed solution, the equilibrium solution (52)-(54) 
will be perturbed by allowing the coefficient of sin Q t to become a slowly varying function of time 
as compared to the period 27ril — 1 of the electromagnetic waves. In addition, it will be assumed that 
ui(p, t) is a separable function (i.e. u\(p, t) = u(p)T(tj). Under this assumption, (85) becomes 

KiSla(p, t) = (f Q (p)+Su(p)T(t))smnt, (Bl) 

where the characteristic frequency of u(p)T(t) is of order w <C fi. This introduces a small time- 
dependent perturbation in the metric functions 

L(p,t) = L (p)+5L 1 (p,t), (B2) 
Q(p,t) = Q (p) + SQi(p,t). (B3) 
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It is stated without proof that for u\(p, t) separable, the field equations impose the relations u(p) = 
fo(p) and T(t) = sin(wt). Therefore the analysis in this appendix is carried out for the perturbation 



Kifla(p, t) = (/o(p) +Sf (p) smut) sin fit, w<fi, S < 1. (B4) 

The perturbation expansion will be taken to the first order in <5. We have assumed that the coefficient 
f(p, t) = fo(p) + Sfo(p) sm °->t of sinfH in (B4) varies on a time-scale much longer than that of 
sinfH. Therefore /(/?, t) is approximately constant over the short time period T = 2irfl~ 1 of the 
electromagnetic waves. Therefore, evaluation of the time-averages (86)-(89) yields 

k? n 2 {a 2 (p, t)) T = \fl (1 + 28smwt) + O (<5 2 ) , (B5) 

K '° 2 ((0) 2 ) =Kd^) 2(1+25sinw<)+0 ^)' (B6) 

K l Q2 ((^) 2 ) =^ 2 .fo 2 (l + 25sin^) + 0(5 2 ), (B7) 



Substitution of (B4)-(B8) into (45)-(48), expanding to first order in 6 and setting each order to 
zero yields the unperturbed equations (63)-(65) and the first order equations (B9)-(B14). Setting 
the first order part of the wave equation (45) to zero yields 

A(p, t) sin fit + B(p, t)cosOt = 0, (B9) 



where 



A( P , t) = n 2 ^ 2 n- 2 /o + (i - i* 2 P - 2 Ql (i - ^)) f + Qo (i ^ ) x 

2o d/o ( T dL \ , ^ d 2 / ^ 2L 



H 2 " ^r-'dfj+^i 1 , 

d/odQo (, 2L \\V,„ , o2 /^ ^ 2L 



+ #T7l 1 -T ; JJj smu ' , + S! A <3 "l 1 -irl> 



x 1 2/9 ~ \ Ql \ L ° ~ + Qo \ Ll ^ t] ~ p d~ P Ll{p > i} 

1 -z-Qi(p, t) - 2p '■L 1 (p, t)- 



p J dp dp J dp 

d 2 ./V 



+ (Qi(p, t) (l - - 2p~ 1 L 1 (p, t)Q 



dp 2 
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+ 2p- 3 ( r 2 QgLi(p, t) - l* 2 P QoQi(p, *) ( 1 - — ) ) /o 



+ (q^, t) (i - ^) - 2p- 1 ii(p, t)g ) (V 2 ^<?o (io - P^ 



2L \ /. d 2 /o , dQod/o 



and 



dp 



B(p, t) = O/o ^p" 1 - ^Li(p,t)-Qo 1 ^Qi(p,t) + 2wcoswt^. (Bll) 



The derivatives of L\(p, t) and Qi(p, i) with respect to p are found from the first order equations 
>o^MP, *) = - ((l - 2 ~y) QiiP, t) - 2p- 1 L 1 (p, t)Q ) ^ 



-ig -g l( p,*)(i(i-^)V„(^) a + ^ 

+ \ (l* 2 p- 2 (l - ^) QiiP, t) 2i*V 3 £i(p, t)Q ) Jl (B12) 



and 



1 - Qo^~ p Qi(p, t) = -(j(l- Qx(p, t) - 2p- 1 L 1 (p, t)Q„) 



dQp 
dp 



p - 2L \\ p ) \dp 



+ \ U (l - ^"ooQ.fo - 4„- (l - 2 -f) £,(,, t)Q*\ 
respectively. The derivative of £i(p, t) with respect to i is given by the first order equation 

s^'J-M'-tO*^ "'* <bm) 
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The simplest of the first order equations is (B14). It is immediately integrable to yield 



Li{p, t) = - I 1 —\ fo— smut + dip), 



(B15) 



where c\{p) is a function of p. It is possible to obtain an equation for dQ\(p, t)/dt from (B9). This 
is done by substituting (B14) into the coefficient of cosfH {viz. B(p, t)) and setting the expression 
to zero. This is justified since sin t and cos t are independent functions in the approximation that 
A(p, t) and B(p, t) are slowly varying functions of time. Solving for dQ\(p, t)/dt and integrating 
yields 



Up to this point the only first order field equation which has been satisfied is (B14). To satisfy 
the first order wave equation (B9), the condition A(p, t) = must be imposed. By using the first 
order field equations (B12) and (B13), the p derivatives of L\(p, t) and Qi(p, t) found in A(p, t) 
(Eq. (BIO)) can be eliminated. After substitution of (B15) and (B16) into (BIO), A(p, t) no longer 
depends on the first order functions L\{p, t) and Qi(p, t). As a result of these substitutions, (BIO) 
takes the form 29 



where A(p) and B(p) depend only on p, the unperturbed functions fo(p), Lq(p), Qo(p) (and their 
derivatives) and the two functions c\{p) and C2(p). Note that C\{p) and C2(p) are only found in B(p). 
Equation (B17) will be satisfied for all t only if A(p) = and B(p) = 0. Since c\{p) and C2(p) are 
yet to be determined, we will focus upon the equation A(p) = 0. 

The function A(p) is comprised of the known functions fo(p), Lq(p) and Qo(p). It is therefore 
necessary to transform A(p) to the x coordinate system and expand in inverse powers of I* 1 / 3 as is 
done for the unperturbed system. 1 The transformation is given by Eqs. (75)— (76). After a lengthy 
computation, the asymptotic expansion of A(p) = yields 



In order for (B18) to be satisfied for large arbitrary I* (in the limit I* oo), each order of Z* 1 / 3 must 
be set to zero. Since setting <p(x) = implies the absence of electromagnetic wave perturbations, 
the bracketed term must be zero. It is known from the unperturbed system that 1 




(BI6) 



A{p)sinwt + B{c 1 {p), c 2 {p), p)=0 



(B17) 




(B18) 



A (x) = -(l-fc 2 (x)). 



(B19) 
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Substitution of (B19) in (B18) leads to the relation 

uj 2 = 4ft 2 (B20) 

which must hold in order for the field equation (B9) to be satisfied. However, Eq. (B20) is a 
contradiction to the original assumption u <C ft. This result is identical to that found for the 
perturbation analysis based upon Eq. (55). 

Appendix C. Angle Average of T£ 

Equations (A11)-(A14) includes integrating over the angle ip and dividing by the solid angle 4n, 
thus all that is left is evaluating the integrals. The dependence of X" M " comes in three forms 

sin" 2 (6,(0)) 2 , sin" 2 9 (0,(0),, ) 2 and sin" 2 9 9,(0)6,(0),,,, (CI) 

where 

Q l (O) = C?B l {0) (C2) 

and 

B,(0) = sin 0-^(^0). (03) 
The exact integrals are evaluated below: 

J\in- 2 9(B l (9)) 2 S in9d9= Hfcil, (04) 

£ sin' 2 9 (5,(0),,) 2 sin0d0 = , (C5) 

P sin- 2 0B, (9)B l (9) M sin A9 = - ^ + 1} . (06) 
The normalization constant for 0,(0) is found by requiring 

/>27T />7T 

/ / |6,(0)| 2 sin0d0d(^ = 1. (07) 
Jo Jo 



Therefore 



2?r \Jo v w/ / 2tt V(2^ — 1)(2Z + 1)(2Z + 3) 

Thus the normalization constant is 



(2Z — 1)(2Z + 1)(2Z + 3) 



1/2 



^ = < v ~" ^v+y 1 " 7 1 • (c9) 
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